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A METHOD FOR STUDYING THE HUNTING OSCILLATIONS OF AN 
AIRPLANE WITH A SIMPLE TYPE OF AUTOMTIC CONTROL 
By Robert T. Jones 



SUIIMARY 



A method is presented for predicting the amplitude 
and frequency, under certain simplifying conditions, of 
the himting oscillations of an aiitomatically controlled 
aircraft with lag in the control system or in the 
response of the aircraft to the controls • If the steer- 
ing device is actuated by a simple right-left type of 
signal, the series of alternating f i:?ced-amplltude signals 
occurring during the hunting may ordinarily be represented 
by a square wave. Formulas are given expressing the 
response to such a variation of signal in terms of the"^ . 
response to a unit signal. A m.ore complex t\npe of hunt- 
ing, which m.ay involve cyclic repetition of signals of 
varying duration, has not been treated and requires 
further analysis. Several examples of application of 
the m.ethod are Included and the results discussed. 



INTRODUCTION 



V/hen an airplane or other aircraft is. directed by 
a simple right-left signal from an automatic steering 
device, the result is usually a maintained hunting oscil- 
lation about the desired path. The. amplitude of this 
oscillation is influenced by the amount of backlash or 
"dead spot^* in the control system, and by the damping of 
the motion of the airplane. In the following analysis 
the a^nplitude and frequency of these oscillations is 
investigated in terms of the response characteristics 
of the airplane. 
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ANALYSIS 



The analy,^!.^ is based on consideration of the 
response of the airplane (in tenns of angle of yaw or 
pitch) to a continued (unit) signal (fig. 1). This 
response may be calculated by the ordinary theory of 
dynamical stability and it will be convenient to repre- 
sent it in operational form (references 1 and 2): 



Rl(t) ^ Ri(D)l(t) (1) 



The unit response ordinarily occurs in the form 



from which is obtained 

r \pt 1 

Rl(t) = C(t) + iCie + ^ + . . . J (2) 

vvhere C(t) is the steady-state motion, C]_ and C2 
are the constant coefficients of the Heaviside expansion, 
and X^, Xg, etc., are the nonzero roots of the char- 
acteristic equation defining the natural periods of 
oscillation and the damping^ of the aircraft without 
signal. The function f (D) and the particular solu- 
tion C(t) depend on the time variation of control 
displacement produced by a signal and on the stability 
characteristics of the airplane in the degrees of 
freedom in which the control operates. (See refer- 
ence 3.) In the case of a continued signal, the usual 
form of the function C(t) is 

C(t) = C^^ + C-ot 

where C^^ is the steady rate of turn called for by the 
signal. (See fig. 1.) 
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During a himting oscillation the automatic steering 

device reverses the signal periodically as the airplane 
swings through the desired heading* A typical hunting 
oscillation is shown in figure 2. Here it is assumed 
that the reversal of signal is delayed either because of 
a "dead spot" in the steering device or because of back- 
lash in the control mechanism or a combination of the 
tv/o. As indicated, the oscillation will have a funda- 
mental period 2Tr/(0 but may also involve components of 
higher frequency, depending on the natural modes of 
oscillation of the airplane* Ordinarily the shorter- 
period com.ponents do not have' sufficient amplitude to 
cause a reversal of the signal during a half cycle. In 
these cases the variation of signal with tim.e will be 
represented by a sim.ple "square wave," which may be 
ej^pressed as a fijyiction of time by 



I ^ n ^^^^^^ ^ 1> 2> ^> ---^ 
n 

or, m.ore conveniently, by the imaginary part of the 
corresponding exponential series; that is, 

n 



Y/here t = 0 is taken to represent a time at which the 
signal becomes positive. 

The response to the alternating signal is obtained 
by substituting expression (4) for the unit function 
l(t) in equation (1)* Thus, 

R(t) = I. P. Ri(D)| e^^^^. • (5) 

n ' 

If the airplane is inherently stable, so that transient 
effects following the start of an oscillation disappear 
with time, the remaining steady oscillation will be 
represented by 

R(t) = I. P. I ^ H Rl(inu))ei^^t 
n 
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Equation (6) give^ the forced oscillation of the 
airplane in response to an alternating signal in the form 
of a square wave cf any frequency 

By investigating the form of these forced oscilla- 
tions at vaiiou?? frequencies it will be possible to 
ascertain v/hether such oscillations, under the condi- 
tions of autcniatic control, will give rise to the 
assumed alternating 3ignalx<5. of equal duration, and thus 
to establish certain ranges of oo over which hunting 
of this type can occur. It will also be possible to 
establish, in these ranges, a correspondence between • 
the frequency of the hunting oscillation and the magni- 
tvide of the dead spot- Vvith the frequency determined, 
it is possible also to find the- amplitude of the oscil- 
lation and the maximutn deviation of the airplane from 
its path* 

In the simplest cases the required information may 
be obtained directly from equation (6) . In the case of 
more complex motions, further analysis will be required 
as follov/s: 

As a first step, separate R-j^(ina)) into its real 
and imaginary, parts 



The functions A and B may be plotted against no) as in 
figure 3. The.«?e functions will show peaks near values 
of"" no) corresponding to the resonant frequencies of the 
airplane. Then, for*^any particular hunting frequency co. 



At the time of reversal of the signal sin noot = 0 
and cos noot = ±1- accordingly as the signal is becoming 
positive or negative. The am.plitude of the response at 
this instant is therefore 



Ri(ina)) = A(na)) + i B(na)) 



R(t) 




cos nw 



a (7) 



n 
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This amplitude will also be the ariplitude of the dead 
spot. (See fig. 2.) A plot of 

Re('^) = I ^ h ^^^^^^ 
n odd 



can readily be obtained from the Evcurve of figure 3 and 
v/111 show the periods of the h\anting. oscillation corre- 
sponding to various widths of dead spot. 

The slope of the response curve at this saine instant 

is 



If the response to a positive signal is negative (as in 
fig. 1), in order that the motion represent a possible 
hunting oscillation (that is, be consistent with the 
assumed variation of signal), it is necessary that 

(I) % I 6 



for a positive dead spot, and that 
(II) R' > 0 



indicating that the airplane crosses the dead spot in 
the proper direction. A further condition is that no 
more than one complete crossing of the dead spot occurs 
v^ithin one-half cycle; that is. 



(III) R(t) > -Rg 

(See fig. 2.) The value of R(t) in the middle of a 
half cycle is relatively simple to obtain;: 

Ra ^ |&(^) - 4 A(3a)) + i A(5co) - • • • J 
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and raay be used as a criterion, though Ra is not neces- 
sarily^^the xriaximura or minimum value of R(t) (see fig. 4) 
and condition III mav not be satisfied even though 

It should be noted that, in the regions excluded by 
the foregoing conditions, a more complex t3rpe of hunting 
oscillation involving a sequence of signals of different 
durations may occur. In these . regions , the R^ and 

curves derived for the square-wave signal no longer apply 
to the condition of automatic control. These oscillations 
require analysis beyond that presented in this report. 



EXAMPLES 



In order to demonstrate and check the procedure 
described, assume a sim.ple response characteristic in 
which the airplane iiiiriiediately starts turning at a - 
constant rate, as directed by the signal. V/ith this 
response 

Rl(D) 
R] (nwi ) 
and, from equation (7), 

17*-" - - 
4 \ 

ri 

which is the Fourier series for a saw-tooth 90^ 
out of pha^e with the signal* (See fig. 5.) In this 
case the response occurs vvdthout lag and the amplitude 
of the hunting is exactly equal to the dead spot. The 
frequency co is ^C|^/2 divided by the width of the dead 
spot. 

A simple exarr.ple nearer the practical case is one 
in which the signal causes a force P to act on a 
mass -m. In this case the response to a \anlt signal is 




—5— cos no:;t 
n"^ w 



and the hunting oscillation is seen to be 



R(t) 




sin nOi>t 



The exriression is recognized as the Fourier series for a 
succession of parabolic segments (fig. -6). It should be 
noted that there is no component out of phase with the 
signal, with the result that Rg is zero for all values 
of CO • Eence the calculation shows no possibility of 
hunting with a finite dead spot. In fact, it can be seen 
from energy considerations that if a dead spot existed 
the oscillation would be divergent. 

Interesting applications of the method are furnished 
b?/ cases in which the response to a signal shows a lag, 
possibly due to backlash in the control mechaniSxTi, in 
addition to a dead spot. A simple example of this l-cind 
is illustrated in figure ?• Kera the response is similar 
to that in the first example (fig. 5) except for the 
time lag r. Use is made of the well-knov/n lag operator 

e-"^. Thus, 



f (t) = f(t - T) 



e 



Aoplying this operator to the response in figure 5 one 
obtains 




\ 
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With the lagging 're.aponae , the himting oscillation 
is not coniirxed to the amplitude of the dead spot and, 
in fact, hunting v^lll occur with no dead spot. It is 
easily seen by reference to figixre 7 that the half 
period of the oscillation in this case (no dead spot) 
is 

Tr/o) ~ 2t. 



Langley Memorial Aeronautical Laboratory, 
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